Drastic fall-off of the thermal conductivity for disordered lattices 
in the limit of weak anharmonic interactions 
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(f*j Abstract 

We study the thermal conductivity, at fixed positive temperature, of a disordered lattice of harmonic 

CN 

oscillators, weakly coupled to each other through anharmonic potentials. The interaction is controlled 
q_) by a small parameter e > 0. We rigorously show, in two slightly different setups, that the conductivity 

has a non-perturbative origin. This means that it decays to zero faster than any polynomial in e as 
e — > 0. It is then argued that this result extends to a disordered chain studied by Dhar and Lebowitz 
[12] . and to a classical spins chain recently investigated by Oganesyan, Pal and Huse 23:. 
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1 Introduction 

It is generally admitted that the thermal properties of solids can be derived from molecular dynamics, 
but it remains to this day a widely open conjecture, at least from a mathematical point of view [9] |13] [19] . 
Indeed, the few hamiltonian systems that can be handled analytically to a large extent, also appear to 
have a very pathological thermal behavior. So is it for the ordered harmonic chain 26J, for the Toda 
lattice (see Section 6.3 in [7]), and for a one-dimensional system of colliding particles [28]. As we will see 
soon, disordered one-dimensional harmonic chains also fall into this category. 

This said, much progress have been recently accomplished in the case of a harmonic crystal that 
is weakly perturbed by means of anharmonic interactions [15] [5] |20j . In particular, a Boltzmann-like 
equation for phonons can be, at least formally, derived in the limit where the interaction vanishes, after 
that time has been properly rescaled. A similar strategy is implemented in |13)|19| to another kind of 
system. In these examples, the conductivity can thus be understood in a perturbative regime. 

Besides anharmonicity, impurities constitute another possible mechanism destroying ballistic trans- 
port of energy. In disordered harmonic crystals, one assumes that the mass of each atom is random 
|27j[llj. The dynamics is integrable and, in dimension one, the conductivity of disordered harmonic 
chains is by now rather well understood |29j [T] . In particular, in presence of everywhere onsite pinning, 
Anderson localization of the eigenmodes forbids any transfer of energy [6] . Though still integrable in the 
mathematical sense, the dynamics is much less understood in higher dimensions. Normal conductivity is 
expected in dimension three if the disorder is small enough |10) . 

We here look at the effect of both disorder and anharmonic interactions on a lattice of pinned har- 
monic oscillators. Three studies have attracted our attention. First, Dhar and Lebowitz [12] conclude 
from numerical experiments that the conductivity of a one-dimensional pinned disordered harmonic chain 
becomes positive as soon as some anharmonicity is added, destroying thus localization. The same conclu- 
sion is then reached by Oganesyan, Pal and Huse [23], which study a chain of classical spins. Moreover, 
they observe a very rapid fall-off of the conductivity as the parameter e > (called J in [53]) controlling 
the interaction goes to zero. Finally, Basko [2] argues that the conductivity in a comparable system can 
be understood as an effect of Arnold diffusion. He also gives an expression for the conductivity, predicting 
a decay with e — > that is faster than any power law (equations (3.6) in [3] with p in place of e). 

Let us describe our contribution. We study a d-dimensional lattice of 1-dimensional harmonic os- 
cillators with random eigenfrequencies, as described in detail in Section [2] For e = 0, they just evolve 
independently from each other, but, for e > 0, each oscillator is coupled to its neighbors through an 
anharmonic potential of strength e. 

Since the rigorous analysis of the conductivity of this system is likely out of reach, we study it in two 
slightly different simplified setups. Let n » 1 be some large integer that we fix. The first simplifying 
procedure consists in following the dynamics only up to very large time-scales of order e~ n , and study 
the conductivity in the strict weak coupling limit e — > 0, as in [5] [2U] [TS] [IH] for n = 2. Alternatively, 
we can perturb the hamiltonian dynamics by a stochastic noise that preserves the energy of individual 



atoms, and possibly models some chaotic behavior of the system. This noise is however chosen such that 
it only becomes sensitive at time-scales of order e~ n . It destroys hypothetical ballistic motion of energy 
only after such long times. These two approaches are very close in spirit. 

Our results are stated in Section p3l Let < (3 < +00 be some given inverse temperature. Let k(e) 
be the Green-Kubo conductivity at inverse temperature (3 (see (p| below) . Theorem [l] deals with the 
first setup. Let m be an integer such that 1 << m < n. It asserts that that e~ m K(e) vanishes in the 
weak coupling limit, no matter how large the integers m < n are taken. We believe this to be a strong 
indication that «(e) = 0(e m ) for every m > 1. Theorem [2] deals with the second setup, and asserts that, 
in that case, k(c) = 0(e n+2 ). It extends, in a simpler model, results obtained in [5] [5], where the noise 
was dominating on the anharmonic interactions. Incidentally, a comparison of our work with [5] and [6] 
reveals that the effect of noise and that of anharmonic potentials can sometimes be very different in some 
respects. 

Since the conductivity is defined at positive temperature, we find it hard to establish a clear connection 
between our work and the numerous existing results dealing with the possible dispersion of a finite energy 
packet at infinite volume (among many other references, let us mention P3][2I][ST] ^ or mathematical 
results, or |22j for numerical ones). Still, in the perturbative regime that we consider, our results actually 
follow from a purely local analysis ; they are based on the lack of resonances in the system, a common 
phenomenon well described for example in Section 2 of |16j . 

Theorems[T]and[2]are proven in Section|4] In Section[5j we show, without providing a full mathematical 
treatment, how to adapt our line of reasoning to obtain similar conclusions for the disordered chain of 
[T2"] and the spin chain of 23 , in agreement with the observations of [53] and the predictions of [3J. 
While the comparison with our model is completely straightforward in the case of the spin chain, if some 
assumption on the distribution of magnetic field is added, at least some technical work would be required 
to get a true mathematical proof for the chain of p~2]. We finally suggest a connection with a chain of 
weakly coupled identical strongly anharmonic oscillators. 

2 Model 

Let N > 3, and let Zjv be the set of integers modulo N. Let also d > 1. We consider a set of one- 
dimensional classical oscillators, with equilibrium positions on the periodic <i-dimensional lattice 1 N - 
The phase space consists thus in the set of points 

(g,p) = (q x ,p x ) xe ^ N g (R 2 ) Nd - 

The Hamiltonian is written as 

H (1,P) = \ Y. {Pl+u 2 x ql)+e J2 U^ + l E E Vfa-Qv) 

x£Z d N x£Z,% x£Z% y-\x-y\i=l 

with the two following definitions. 



First, the eigenfrequencics of the decoupled harmonic oscillators, uj = (uj x ) xeIi d , form a sequence of 
independent and identically distributed random variables, with law P* independent of N. Let E*(-) be the 
associated expectation. While more general distributions could likely be considered, we will assume that 
P* admits a bounded density, so that in particular E^cj 2 .) — E*(wj;) 2 > 0, and that, for some constants 
< cj_ < u + < +00, the bound cj_ < u x < ^>+ holds almost surely. 

Second, the coupling constant e is strictly positive. The potentials U and V belong to C£f mp (R), the 
space of infinitely diffcrcntiablc functions with polynomial growth. The pinning potential U is taken such 
that, for every a > 0, for some constant c > 0, and for all x G R, 



e 



- a ("-v 2 +eU(y)) dy < +OQ and oj_+eU"(x) > c. 



The interaction potential V is supposed to be symmetric, meaning that V(—x) — V(x) for all x € R. 
Moreover, in order to ensure good decay properties of the Gibbs measure, we require the pinning to 
dominate over the interaction: similarly to (2.3) in [8], one asks that there exists a constant C < +00 
such that, for every x,y £ R, 

V"{x-yf < C(l + U"(x) 2 )(l + U"(y) 2 ). 

We possibly perturb the hamiltonian dynamics by a stochastic noise that preserves the energy of each 
oscillator. Let e' > 0. The full generator of the dynamics writes 

(0) , ,4(1) 



with 



L = A har + eA anh + e'S = A har + eA ( ^ h + eA ( J h + e' S 

AharU = ^ (Px 9q x U - wlq x d P:c u) , 

41« = -E u '^) d v^ 

x&% 

Su = ^2 (u(---,-Px,---)-u(...,p x ,...)). 



We denote by (X^ ) t>0 , or simply by (X t ) t > , the Markov process on (R 2 )^ generated by L. 

Besides E*Q, we will consider two other expectations. Let < /3 < +00 be an inverse temperature, 
that will be considered as a fixed parameter in the sequel. Let {■) p be the corresponding Gibbs measure, 
defined by |2| below. We will just denote by E(-) the expectation over the realizations of the noise. We 
will write E / g(-) for (E(-)),g. 

We need some extra informations on the Gibbs measure. Its density pp is given by 

P p(q,p) = e-W**)/Z(J3), 



where Z({$) is a normalization constant. This density factorizes as pp(q,p) — f/g{q) . p'Lip), and the 
density p'a factorizes again: 

m - (lf n n *-*'»■ 



xQ 



N 



While p'g(q) does not factorize, the conclusion of Theorem 3.2 in [H] holds thanks to our hypotheses on the 
hamiltonian .HI 1 ! there exist constants C < +00 and c > such that, for any functions u, v £ C^ mp {R N ) 
satisfying (u)p = {v)p = 0, it holds that 

\if.g) fi \ < Ce-^(«)^))((v gU ) 2 )y 2 ((v 9 ,) 2 )y 2 . (1) 

Here, S(u) is the support of a function u, defined as the smallest subset of 1 d N such that u can be written 
as a function of the variables q x for x £ S(u), whereas d(S(u), S(u)) is the smallest distance between any 
point in S(u) and any point in S(v). 



3 Conductivity 

Let x £ Z^r . The energy of particle x is given by 

e x = -{pl + v x q 2 x )+eU{q x ) + - ^ V (<lx~qy)- 

y:\x—y\i=l 

For y £ Zfj- such that \x — y\\ = 1, one defines the current 

Px+Vyrri, V 

ej^y = e g ^(& _ 9»)- 

Since V is symmetric, one has j y<x = —jx.y With this definition, one computes that 

d 



Lie x e J \Jx—e k ,x ]x,x+e k/ 



fe=l 
Uth 



where e^ is the k unit vector. 

To simplify some further notations, we find it convenient to define now a set of oriented bonds in the 
lattice 7L d N : 

Z d N = {(x,y) £ (7L d N f : y = x + e k for some 1 < k < d}. 

The total current and the rescaled total current are defined by 

eJ N = e ^ Jx and eJ N = -ttjj^Jn- (2) 

x&% 

Then, if the limits exist, the Green-Kubo conductivity of the chain is defined by (see for example [IB] ) 
K(e,e') = lim lim «(e, e', N, t) = lim lim d -1 ^ 2 Esf 4= / Jn ° X^' €>) ds ) . (3) 

t— s-oo AT— j-oo t->oo iV->oo \\t Jo J 

1 Strictly speaking, a slight adaptation of this theorem is needed to cover our needs: going through the proof, one checks 
that the randomness on the pinning does not affect the result, and that the constants can be taken independent of e. 



It maybe could look more natural to define the total current, and the associated Green-Kubo conductivity, 
with respect to a given direction, say e^, for some 1 < k < d. This means that the sum in (pi) should be 
restricted to currents on bonds of the type x — (x, x + e^) for x G Z^. Our results below still hold with 
such a restriction. 

We first prove a result on the fall-off of the conductivity in a weak coupling limit ; it is physically 
most meaningful for 1 << m « n: 



Theorem 1. Let n > 1 and let e' = 0. Let also 1 < m < n. For almost all realizations of the 
eigenfrequencies, it holds that 

imsuplimsupe~ m ( [ — ^= / 

t— >oo 



lim lim sup lim sup e^" 1 / [ —L= f J N o X^ ds ) \ =0. 



If the energy conserving noise is added, the following also holds: 

Theorem 2. Let n > 1 and assume that 

e' = e n . 

Then there exist constants C = C(n) < +oo and e > such that, for almost all realizations of the 
eigenfrequencies, and for all e £ [0, eo], one has 

ft \ 2 



lim sup lim sup E 



p(A=J JnoX^As^ < C(n).e^ 



Remarks. 1. The hypothesis that the law of the frequencies (u> x ) xeZ d admits a density cannot be 
completely dropped out for the proof to work. Assume indeed, at the opposite, that the frequencies 
should take only two values. In that case, among three oscillators, two at least have the same frequency. 
Then, as long as one takes d = 1 and n > 1, or d > 2, some fatal resonances appear and invalidate our 
line of reasoning. 

2. It is possible that the weak coupling limit or the noise actually destroys some effects, such as a very 
slow ballistic transport, that become dominant when looking to the system at large enough time-scales. 
For example, let us assume that n is given, that d — 1, that all the interactions are harmonic, and that 
(uJ x )xez N form a sequence of 2n different frequencies, arranged periodically. One checks that our proof 
works in this particular case, although the true conductivity is known to be infinite. Since however, in 
our theorems, n can be taken arbitrarily large, such a scenario only becomes possible at time-scales larger 
than any inverse polynomial, and does not affect our main message that the transport of energy has a 
non-perturbative origin. 

3. The proofs are achieved by gathering three different ideas. Given x € Z^ , we first approximately 
solve the Poisson equation — Lu = j x , assuming a non-resonance condition on the random frequencies 
of the oscillators located near x, as defined in the beginning of the next Section. We next exploit the 
conservation of energy to see that the currents j x , at places x such that this condition is violated, in fact 
do not contribute to the evaluation of the conductivity. We finally take into account that e — > before 
that t — > oo for Theorem [T] or that S generates a diffusion on time scales of order e™ for Theorem [2] in 
order to get a bound on rest terms. 



4 Proof of Theorems [T] and [2] 

The proofs of our two results are analogous but, since that of Theorem [2] is slightly more involved, we 
focus on it. The needed adaptations to show Theorem [T] are presented in the last Subsection. 

Let us first define what are non-resonant frequencies near a bond x G Z^. Let £ = (£, x ) x ei, d w hn 
£r G Z for all x G Z^. Given an integer r > 0, given x G Z^ and given x = (a, 6) G Z^, we define 

Bi(s,r) - {yGZ^: ||/ — as|i < r} and B x («,r) = B^o.r) UB x (6,r). 

Given x G Zjy-, given a realization w of the frequencies, and given an integer r > 0, we define 



j/GBi(x,r) 



For a, c G R, for x G Z 1 ^, and for an integer r > 0, we next define the diophantine set 

D a , c (x,r) = ju; : |(w,0£,r| > jTT^ for all ^0 on Bi(x, r)|. 

Here, ^/OonB^^r) means that £ x ^ for at least one x G Bi(x, r). The following monotony property 
holds: if r < r', then D atC (x,r) C D a>c {x,r'). It is known (see for example [23]) that, given 8 > and 
an integer r > 0, there exist constants a < +oo large enough and c > small enough such that, for all 
x G Z^, , it holds that 

P*(weA,, c (x,r)) > 1-6. (4) 

From now, let n > 1 be as in the hypotheses. Let S > to be fixed in Subsection [42] and let a, c > 
be such that Q holds with r = n. We also fix a realization u of the frequencies. All the constants 
introduced below may depend on n, <5, a and c. 

4.1 Approximate solution to the Poisson equation 

Let us start by a simple, but crucial, observation. We say that a function u — u(q,p) on the phase 
space is p-symmetric if u(q,p) — u(q,—p) for all (q,p) G (R 2 )^ ; we similarly define p-antisymmetric 
functions. The generators Ahar and A an h map p-symmetric functions to p-antisymmetric functions, 
and p-antisymmctric functions to p-symmetric functions. The interplay between p-symmetric and p- 
antisymmetric functions has shown to have deep consequences in an other context [J]. 

-i TTk/Tn,2N d , 

[■)/3> 



For the rest of this Subsection, we fix x G Z^, and we assume that ui G D ac (x, n). For an integer 
k > 0, we designate by ||| • |||fe the norm of the Sobolev space H fe (E 2A ' , (•) / g) ; explicitly 



2iV" 



2N d 

INIU = Y, W d( * u \\v>(m™ d , {■}?,) with \ a \i = J2 a J and dau = d xl---df 2Nd - 

aeN 2 « d :|Q| 1 <fe 3 = 1 

Let us first show the two following lemmas. 

Lemma 1. Let ui G D ac (x,n). Let f G C£f (R 2Ar ) be a p- antisymmetric function that depends only on 
the variables (q y ,p y ) with y G Bi(x, r) for some < r < n — 1. Then there exists a p-symmetric function 



U G Cj^ mp (R 2N ), depending on the same variables as f, satisfying (u) p — 0, and solving 

~Ah ar u = j. 
Moreover, there exists a integer m > such that, for any k > and for some constant C k < +oo, 

HHIU < c fe |||/||| fc+m . 



It follows from this Lemma that A an -^u € C^ (M. 2N ) is p-antisymmetric, and depends only on the 
variables (q y ,p y ) with y £ Bi(5;, r + 1). Moreover, there exists some non-negative integer m' such that, 
for any k > and for some constant C k < +oo, 

|||Awi/iw|||fc < Cfe |||/|||fe+ m '. 

Proof. For z > 0, the solution u z £ C^ mp (R 2Nd ) to the equation 

(z - A har )u z = f 

exists and is unique ; it is given by 

/>oo />oo 

u z = e- zt e Ah ^fdt = / e~ zt foX? ar dt (5) 

Jo Jo 

where (X' t lar ) t >o is the deterministic process generated by Ah ar - Since / £ C^ mp (R 2N ), and since 
{X^ ar ) t >Q explicitly given by ([7]), it is checked that u z defined by (|5| is indeed an element of C^ mp (R 2N ). 
We will show that there exists a function u, having all the regularity properties appearing in the conclu- 
sions of the Lemma, depending on the same variables as /, and such that 

lim u z = u and lim A har u z = A har u, (6) 

0<z->0 0<z^0 

these limits being for example understood in L 2 (E 2JV ). This implies that u solves the equation —Ah ar u — 
f . So does then the function u — (u)fj, which inherits also from the other properties of u mentioned in 
the conclusions of the Lemma. This will thus conclude the proof. 
For given initial conditions (q,p) G M. 2N , one has 

x t ar {^P) = (• ■■,qkC0suj k t+ (p k /uj k )smuj k t,-uj k q k smuj k t + p k cosuj k t,. ..). (7) 

Therefore, defining 

g(q,P;9) = f(---,Qk cos Ok + (j>k/uk) sin 6k,-0Jkqk sin Ok+Pk cos Ok,---) 

and 

0(t) = LoteR Nd , 

one has 

/>oo 

U z (q,p) = / e- zt g(q,p;0(t))dt. 
Jo 



It is convenient to write g in Fourier's variables: 

g(q,p;8) = V gfapiQe* 6 * with $(g,p;£) - — \^ f g{q,p; 9) e~ i9 < dO. 

Since / only depends on the variables (q y ,p y ) with y £ B(x, r), it holds that g(q,p,£) = for any (q,p), 
as soon as £ y ^ for some y ^ B(x, r). Let us see that g(q,p; 0) = for any (q,p). For this, let us define 

9 = 9(q,p)e[0,27T] Nd by 

q k +ip k /u k = (ql+pl/ujl) exp(i9 k ). 

The p-antisymmctry of / implies that 

g(q, P ;29-9) = -g(q,p;6), 

from which g(q,p; 0) = follows. 
One now computes that 



/>oo 

and so, taking the limit z — > 0, 



g(g»p;0 

z — iui ■ £ 
£ez« d -{o} " eez« d -{o} 



«(«)P) = » zl 



g(q,p;0 

£SZ« d — {0} 

The regularity properties of u are obtained thanks to the hypothesis ui £ D ac (x,n), as is usual in KAM 
theory ; see for example |25j . From there, ffity follows. D 

Lemma 2. Let f £ C^ m „(R ) be a p- antisymmetric function that depends only on the variables 
(q y ,p y ) with y £ Bi(x, r) for some < r < n — 1. TTien i/iere exists a function u £ C£? mp (R ), 
satisfying (u) p = 0, depending on the same variables as f , and solving 

-Su = f. 

Moreover, for any k > and for some constant C k < +oo, it Zio/ds that \\\u\\\k < Cfc |||/|||fc. 

Proof. Let us fix the values of |p x | and g^ for all x £ 1%, so that in fact / is now seen as a function on 
{ — 1, -f l}l- B i( :E > , ')l , the set of signs of the impulsions. The operator S generates a standard continuous time 
random walk on this set. Writing S = \B\(x, r)\(T — Id), one has u = , B r~ r)\ Sfe>o T k f . The terms T k f 
converge exponentially fast to 0, provided that / is of mean zero with respect to the uniform measure on 
Bi(x, r). So it is since / is p-antisymmetric. The regularity of u follows form this representation. □ 

The function j% is p-antisymmetric, belongs to C^ mp (M. 2N ), and depends only on the variables in 
B(x, 0). Lemma 1 ensures that there exist functions u~ , . . . , u~ solving the following finite hierarchy 

-A har u { ~ ] = jx, (8) 

-A^uf = A anh U { }\ (9) 

-A har uf = A anh 4 l - 1] . (10) 



By Lemma 2 there exists also a function i>~ such that 

OV~ — AanhU- ■ 



Defining 



one finds that 



4 1 ' + eu¥ ] + ■■■ + e"- 1 ^ and v 7 = u 7 - 4" } 



» 



j s = -Lux + e n (A anh u { ~> +Sux) = -Lux + e n Svx. (11) 

The functions u^ and u^ belong to C°°(M. 2N ), are of mean zero, depend only on the variables (q y ,p y ) 
with y £ Bi(x, n), and have (Sobolev) norms that do not depend on x, as long as u> £ D aiC (x,n). 
Remark. While it is not difficult to solve the finite hierarchy ([8])-(|9]) as we did under a non-resonant 
assumption on u, one could in fact have expected it to be much less straightforward. Indeed, the dynamics 
generated by Ahar is very degenerated since it preserves the energy of every single particle. This has two 
consequences. First, the equation —A^arU = / can only be solved if / is of mean zero with respect to every 
microcanonical surface defined by the value of the energies of each uncoupled particle. Second, if this is 
the case, then u is not unique, since every function depending only on the energies of decoupled particles 
lies in the kernel of Ahar- The fact that the currents are p-antisymmetric, and that A an h interchanges p- 
symmetric and p-antisymmetric functions, guarantees therefore that the zero-mean condition is satisfied 
at all orders, without that we have to select a particular solution at each step. The operator S instead 
destroys this symmetry of the system, allowing diffusion of energy. 

4.2 Resonances 

We fix a realization to of the frequencies. We let 

A N = {xeZ d N ■. oj e D a , c (x,n)}. 

We then partition 7L d N — An as 

T^n ~ An = [JC 7 , 

7 

where (C 7 ) 7 is the set of connected components, or clusters, of 7L d N — A^ . We let i(C 7 ) > 1 be the length 
of the longest self- avoiding path in C 7 . Given x £ Z N — An, we also define C(x) as the cluster C 7 such 
that ieC r 

In the sequel, by a slight abuse of notation, we will say that x £ Z N belongs to E C Z N , if there 
exists a bond x £ E such that x is a component of x. So any subset E c 1 N wu ^ a ^ so ^ e considered as a 
subset of Z N . Given K. C Z%, we define its border dJC C Z N as as the set of bonds x — (a,b) such that 

dfC = {x= (a,b) £Z d N : (a £ JC and b<£ K) or (a g K, and b € £)}. 

Given K, C % N) one defines dK, as the border of JC seen as a subset of "Z N . The following Lemma contains 
the core of the argument that allows us to get rid of resonances. 

10 



Lemma 3. Let C 1 C Z^, — An be a cluster, and let £(C 7 ) — m. It holds that 

E •?* = E 0i(^>7)ix - e_1 ^( E ^(x)e x 

xec-y xedC-, xeCyCi.% 

where 9\ and 02 are integer valued functions that satisfy the bounds 

sup \O x {x,i)\ < Gm d , sup 1 02 (a;) | < Cm d , 

for a constant C that depends only on the dimension d. 
Proof. Given K, C 1% , the conservation of energy implies 

e-^C^e.) = ]T Q{x)k. (12) 

Here, writing x = (a, b), one has C,{x) = 1 if b € K and (,(%) = —1 if a G /C. To apply this formula, we 
can repeatedly consider contours that partition C 7 into two pieces, and get rid of all the currents in C 7 . 
A systematical way to do this is as follows (see Figure [I]). Let us give a direction k G {1, . . . , d} and a 
value s G Zjy such that for some x = (a,a + e k ) G C 7 , it holds that a k — s, where a k denotes the k th 
coordinate of a. We then define 

P(s, k) — {x = (a, a + e fe ) G C 7 : a k = a}, 



as well as the set IC(s, k) c C 7 C Z^ made of points x such that x k > s. Then, using ( 12 1, 



e lL [ E e 'A = E ix + E Cfflte- 

^xeK(s,k) ' xeP(s,k) xedK(s,k)-P(s,k) 

It hols that dlC(s, k) — P(s, k) c dC 1 . This operation can be repeated with every direction fc G {1, . . . , d} 
and every value s G Zjy. This only needs to be done 0(m d ) times to cover the whole cluster ; summing 
one obtains 

-e _1 -M E e 2(x)e x j = ^2 jx ~ 51 0i(*>7).?*) 
^xec 7 ' xeC 7 zedc T 

for some functions 0i and 02 that satisfy the stated bounds. □ 

Thanks to this Lemma, and using the fact that L(e)p = 0, we decompose the total current as follows 

Jn = Yrx + Y(Y^) - E rx + E ( E *(2)*) - ^(E ( E «»(*)«.)) 

xSAjv 7 xGC T x£A n 7 a;S9C T ^ 7 x£C-, ' 

=: ^ Vi(^)i2 - e-'lf £ Mx) (ex - (e)p)\ (13) 

x£A N xeif, 

Here 

^i(2) = 1+ E *i(s.7). 

the sum containing two non-zero terms at most, and 

■02 (%) — 02 (^ if xgI JC 7 and ip2{x) = otherwise. 
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Figure 1: The sets C 7 , P(s,k) and IC(s,k) appearing in the proof of Lemmal3J Here d = 2 and 
k — 1. The whole figure represent a cluster C 7 . For a given s, the bonds in red constitute the set 
P(s, k = 1), while the bonds in blue constitute the set IC(s, k — 1). 

The functions ipi an d ^>i inherit the properties of Q\ and # 2 stated in Lemma pi If x ^ V)~dC 1: then 
|^i(a?)| = 1, whereas if x £ {J~dCy, then |^i(x)| < Cm d , where m = f(C 7 ), where C 7 is the largest 
cluster such that x £ 9C 7 . Similarly, if x (jz. IJ-,^, then |^2(x)| — 0, while if x £ C 7 for some 7, then 

hMaOl < Cmd > with m = e ( C i)- 

We now may fix the parameter <5 introduced at the beginning of this Section. We take S > small 
enough so that there exist constants C < +00 and d > such that, for every integer m > 1, and for 

every x £ Zfj, it holds 

P*(e(C(x))=m) < Ce- C ' m . (14) 

The number 5 is well defined. To see this, we first notice that, for x,y £ Z^, the events uj ^ D a ^ c (x, n) 
and w (fi D atC (y,n) become independent as soon as B\(x, n) n Bi(y, n) = 0. Next, we observe that the 
number of self-avoiding paths in 1 d N of length m is bounded by (2d — 1)™\ But the number of disjoints 
balls of fixed radius n in Zfj, with center contained in any such path, is proportional to m as m — > 00. 
Therefore, the probability of forming a path of length in in 1 d N — Am decays exponentially with 771 if S is 
small enough. 



Lemma 4. For every x £ Z N , we let f% be a smooth function on the phase space M. 2N , that satisfies 
{fx)p = 0, that depends only on the variables in B(ir, n), and which iV^-norm can be bounded by a constant 
independent ofx. Then, there exists a constant C < +cxd, such that, almost surely, 

2\ 

- \ ' ,/,, (Z\ f~ ' 

/ 1 



x£A N 



< c, 



lim sup 

JV->oo 



\N d / 2 



]T M*) (^ - (e)*) ) < C 



(15) 
(16) 



Proof. Throughout the proof of this lemma, C < +00 and d > denote constants. The proof of both 
bounds are analogous, and we only handle the first one. Expanding the square and applying Cauchy- 
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Schwarz inequality yields 

Vn := ((~ E Mx)h) 2 ) < ^ E ^i(^i(y)KM^I < ^ E ^?(2)K^/»>/»l 

Then, first, thanks to the bounds on |-0i(x)| stated after the proof of Lemma [31 and thanks to the 



exponential bound ( 14 ) on the probability of forming clusters of size to, it holds that 

E*(|Vi(5Qr)<C p , (17) 

for any p > 1, and for some constant C p < +oo independent of x € Z^r. Second, thanks to the 
decorrelation bound (IT]) on the Gibbs measure, and to the hypotheses on the functions fy, it holds that 

P* := E \(kfv)p\ < C - 
y£A N 



Therefore 



Vn ^ Wd E W(2)-E*W^)))«?+ 1 J S E E *(^(^))/% - >V + ^. 



Since Zm is a bounded deterministic sequence, it suffices now to bound Yjv- By the Borel-Cantelli 
lemma and Markov inequality, the proof will be concluded if one finds p > 1 such that 

E E*(l^| p ) < +oo. (18) 

JV>1 

If d = 1, one can take p = 4, but if d > 2, the choice p = 2 suffices. We only will deal with this second 
case ; the case <i = 1 is obtained similarly at the price of slightly longer computations. To simplify some 
further notations, let us write 

4>(x) = (Vftaf)-E,(^(2))). 

One has 

E*(3#) < ^ E E*(^(i)^(5)). (19) 



For x,y £ 1%-, let us set |x — t/|i = min{|x — y\\ : x £ x,y £ y}. By (17), the terms in the sum 



in the right hand side of ( 19 1 are uniformly bounded, and it now suffices to establish decay bounds on 
off-diagonal terms. Let us fix some x, y € Ajq. The claims below are valid for \x — y\\ large enough. 
Given c" > 0, we then define an event E as follows: u> € E if nor x nor y do belong to the border of a 
cluster of size larger or equal to c"\x — y\. Then, if c" small enough, the variables 4>{x) and <fi(y) become 
independent with respect to P*(-|i£). Therefore 

E*(</>(x)<t>(y)) = E*(<t>(x)<f>(y)\E)P*(E) + E^x^^E^P^E^ 

= E*(^(£)|B)E*(^(jf)|^)P*(^) + E,(0(5?)0(y)|i? c )P,(ii; c ) 
= -E*(^|£)E*(<Ky)|£ c )P*(£ c ) + E*(0(£)<K£)|£ c )P*(£ c ), 
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where one has used the fact that E*(0(y)) = to get the last line. The exponential bound P*(-E c ) < 
Ce _c l^"^ 1 holds, so that in particular P*(-E) is bounded form below by a strictly positive constant. So 
it holds that 

\E*(<t>(x)\E)\ < C|E*(0(S))| < C, 

\E^(y)\E c )P+(E c )\ < £^ 2 {y)) 1/2 Pl /2 {E c ) < Ce" '^ 1 , 
\E4^Wy)\EC)P+(E<)\ < E^ 2 (i)</> 2 (y)) 1/2 py 2 (£ c ) < Ce' c '^'^. 

Thus 

|E*(0(£)<Hy))| < Ce-^l^l 1 . 



Inserting this estimate in (19), it follows that E*(YJy) < C/N d . Since we have assumed d > 2, this in 



turn yields (18), and so concludes the proof. □ 



4.3 Concluding the proof of Theorem [2] 

We fix a realization lo of frequencies. We combine (jTTI) and (fl3| to get 



Integrating over time and dividing by t x l 2 N d l 2 , one finds two martingales M t and 7V t such that 

-^ y ^ o x s ds = m 4 + fl/2Nd/2 E ^) ( W£ - u - ° x *) 

e -l 

+ N t + t i/2 N d/2 E Ms) (( e * - (e)/») ° *t - (ex - (e)/?)) 






xe 



/ o 5 (]vd72 £V^(«)«3f ^.jd«. 



We now separately estimate the variance at equilibrium of the five terms in the right hand side of this 
equation. We remind that, by the result of Subsection |4.1| the zero-mean functions u^ and uj are smooth 
and local, uniformly in x € An- All the estimates below hold almost surely in the limit N — > oo. 
Let us start with the two martingales. First, by Il5|), 



E ^( M * 2 ) =((j^J2 ^(2)«s)(^_E MS)* 



x£A N xEAjs 

o i in *? 1 /O 



Next, since Se^ = for every x G Z^, 



E ^K) = ({mi E ^(x)(e x - (e) fi )\ (-^ E *W(^-(^))) = 0. 



Let us then consider the two terms not involving the operator S. By ( 15 ) and ( 16 ) and the stationarity 
of the measure (■) p, 

V xGAjv / * V xeA N ' ' P 
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and similarly 












E*( 


,-i 


e x o X t ] 


J = 0(c" 


H~ l ). 




L/ V 


/2 N d/2 z^^ 2 ^ \ Cx 

x&% 




One finally deals with the last 


term. By a classical bound [T7] 


there exists 


a universal constant C < +00 


such that 














dsj 








* a2 i( N % 


x£A N ' 


if s 


- 51 V>iOh 

x£A N 


')>, ■ 


= 0(e n ), 



where (15) has been used to get the last bound. The proof is concluded by taking the limit t — > 00. □ 



4.4 Proving Theorem [T] 

Since the proof of Theorem IT] is much similar to that of Theorem [2l we only indicate the main steps. 
We let n > 1. It is enough to consider the case m = n ; we then need to take the limit of the following 
quantity: 

As in the proof of Theorem [TJ we let 6 > as fixed in Subsection |4.2| and we let a,c > be such that 
Q holds with r = n. We also fix a realization uj of the frequencies. Assuming that u> e D a>c (x,n), we 
write as in Subsection |4.1| that 

jx = -Liix + e n A anh uy =: -Lu^ + e n w x . 



All the results of Subsection 4.2 still hold, and we combine this with the decomposition ( 13 ) of the total 
current, to get 

Til jNds= 7ll wr*J^^ u * ds + 7iL m* ^ >(*)(«- -^^ 

x£A N 16ZI 



£ n+l /•£ "t 



/ #5/2 E V'lW^ds. 

x£A N 



We then obtain bounds valid almost surely in the limit N — > 00. Since the generator L is antisymmetric, 
the variances of the terms involving the generator L read 

I wkJ2 M*)ux) ) = ((^j^_£ Mx)(u~ x oX e -, 

= 0(e 2 t- x ), 

r e~ n t 



t,L ^S Wi)e -) >. - or,) 



a:^A N 
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Vf 



Finally, by Jensen's inequality 

= et ((wN E MQwm))^ =0(et). 

The proof is achieved by performing successively the limits e — > and then t — > oo. D 



5 Related models 

We show how to extend our result to two different physical systems, and suggest an analogy with a chain 
of anharmonic oscillators. We only deal with the weak coupling limit, and clearly indicate where some 
pieces of information are missing to get true mathematical statements. 

5.1 Weakly perturbed one-dimensional disordered harmonic chain 

In |12j . a one-dimensional disordered harmonic chain perturbed by small anharmonic interactions is 
considered. The phase space is the set of points (q,p) € M. 2N , and the hamiltonian is of the form 

H &p) =lJ2(pl+ <■&£) + \ E («» - <^) 2 + \ E it + \ E («* - <^) 4 ( 2 °) 

x^TLn a;GZjv zGZ^v xG^n 

= H har + (H anh . 



The hypotheses on (wj;)^ are as in Section 2r The time evolution is predicted by the usual hamiltonian 



equations. At e = 0, the chain is harmonic, and any state can be written as a linear combinations of the 
eigenmodes (£ )i<fc<jv of the systems ; moreover, due to everywhere onsite pinning and since the chain is 
one-dimensional, all these modes are typically exponentially localized. See for example [5j and references 
therein. Let us also denote by (tJk)i<k<N the eigenfrequencies of the modes, which, as opposed to (ui x ) x , 
do not form a sequence of independent and identically distributed variables. At least in a naive picture 
of Anderson localization, it should be however reasonable to think that the frequency wi of a mode £i 
localized near a point x\ G Zjv, and a frequency uji of a mode £2 localized near a point xi € Zjy, quickly 
decorrelate as \x\ — xi\ grows. This should be needed to get rid of resonances as we did in Subsection 



In this model as well, the generators Ahar and A an h exchange p-symmetric and p-antisymmetric 
functions. Now, the crux of the matter is that, for a function / that is p-antisymmetric, and that only 
depends on variables around some site x, the Poisson equation 

— Ahar u = J 

should be solved in the same way as we did, with a solution u that is exponentially localized near x. Let 
us mention that this equation has been solved in (Bj in the case where / is the current corresponding 



In |12| , randomness is on the masses and not on the pinning. This however should not make any crucial difference. 
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to a harmonic interaction. This case is quite special however, since there, a lucky cancellation makes 
disappear all possible resonances. 

Looking at the proof of Lemma [I] one has indeed 

/>oo 

u(q,p) = lim / e- zt foX? ar (q,p)dt. 

This can be made more explicit since, using (4.4) and (4.5) in [5], one can write 

/ N 

foX? ar (q,p) =f ...,^(( g ,e>cos^< + ((p,^)M-)sin^i), 

Here as well, one can define g(q,p; 9) by replacing ut by in this last formula. Defining then an angle 9 
such that 9j is the phase of (q, f ) + i(p,^ J }/uJj 1 one then finds that the p-antisymmetry of / translates 
into the relation g(q,p,29 — 9) = —g{q 1 p,9), which ensures that the Fourier expansion has no constant 
mode. Finally, the exponential localization of the eigenmodes should guarantee that only a small number 
of eigenfrequencies are practically involved, so that diophantine estimates can be used. Moreover, while 
the solution u itself will not be local anymore, it will be exponentially localized, and the limit e — > 
should allow us to get rid of exponential tails. 



Remark. For the chain defined by (20), it is known [2 that the Green-Kubo conductivity n(/3,e), seen 



as function of the inverse temperature f3 and the coupling strength e, satisfies the exact scaling 

k(/3, e) = K,(/3/r, e/r) for any r > 0. 

Therefore, the behavior of the conductivity in the small coupling regime corresponds to its behavior at 
low temperatures. Many numerical works deal with the diffusion of an initially localized energy packet 
(diffusion at "zero temperature"). The heat equation can sometimes serve as a good phenomenologi- 
cal model to describe the spreading of this packet [35]. If one assumes that the diffusion constant in 
this equation is given by the Green-Kubo conductivit}! 3 ! then our results predict in fact a subdiffusive 
spreading of energy slower than any power law. 

5.2 Disordered classical spin chain 

In [53], a one-dimensional chain of classical spins is studied. Let S be the unit sphere in R 3 , and let a 
point of the phase space be given by S = (S x ) x ^z N G § W - The hamiltonian H is given by 

H(&) — 2—i ^k • S x + e 2_^ S x ■ Sx+i, 

and the equations of motions read 

S x = V X HAS X , (21) 



3 While this looks reasonable, it is fair to say that this identification is not justified, since we do not look at the system 
in a true macroscopic scale. 
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where V^ = (d/dSi tX ,d/dS2, x ,d/dS3 >x ), where S • T is the standard scalar product in R 3 , and where 
S A T is the standard vector product in K 3 . It is assumed that (u> x )x£Z. N form a sequence of independent 
and identically distributed random vectors in IR 3 . 



The dynamics (21 1 is generated by A, defined as 



Au = J2 (V X #AS X )-V x u = J2 ( w « A Sx)'V«« + e £} ((S x _i + S x+ i) A S x ) • V x u 

xgZjv iGZjv KgZ N 

= -*4-/iar u + e A an hU 

where the denominations Ahar and A an h are aimed to bare the analogy with our model. The energy of 
the particle at site x is defined as 

e x = u x • S x + -S x _r • S x + -S* • S x+ i, 

and the relation 

Ae x = e(i x _i, x - j x ,x+i) 

is satisfied if one defines the current j x , x +i by 

jx,x+i — 7. • (S x A S x -|_i) + ~((S x _i A S x ) • S x +i + (S x A S x _|_i) • S x +2j- 

At e = 0, each spin just precesses around the axis characterized by the vector ui Xl at a constant 
angular velocity |tt> x |2. Its energy uj x ■ S x is conserved, and determines the plane perpendicular to u> x 
where the precession takes place. The picture in phase space is thus here completely analogous to that 
of the one-dimensional harmonic oscillators we have considered. 

To proceed as before, we would now need to partition the set of functions on the phase space into two 
subspaces, playing the role of p-symmetric and p-antisymmetric functions. This is easily made possible 
if we add the extra assumption that, writing oj x = {uJ\ tX ,'jJ2,x,^>3,x)i one has almost surely 

uj^ x = for all xeZ N . (22) 

We then define 1-symmetric functions on S w as functions that are symmetric under the reflection that 
simultaneously maps io\ iX to — uj\^ x for every x £ Z^v ; 1-antisymmetric functions are defined similarly. 
Then the currents j x ,x+i are 1-antisymmetric functions, and the generator A maps 1-symmetric functions 
to 1-antisymmetric functions and vice versa. This is all what is needed to compute a finite number of 



perturbative steps as we did. Without the assumption (22 1, we have not been able to find a symmetry 
that guarantees the solvability of the hierarchy ([8])-((9]). Still, such a symmetry is noway a necessary 
condition ; we believe that this is only a technical question, and that these equations can in fact be 
solved. 
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5.3 One dimensional chain of strongly anharmonic oscillators 

Wc finally suggest an analogy with a a one-dimensional chain of identical strongly anharmonic oscillators. 
The phase space is the set of points (q,p) £ M. 2N , and the hamiltonian is given by 

Equations of motion are the usual Hamilton equations. For such a system, existence of breathers is well 
known in the infinite volume limit [21j . Unfortunately, this as such does not say us much about the 
conductivity of the chain. 

Let us see that, at least in the asymptotic regime we consider, the effect of the strongly anharmonic 
onsite potentials could be compared to that of random frequencies. At e = 0, all the oscillators are 
characterized by an energy e x .o, and evolve independently from each other. Since their dynamics is 
one-dimensional, they are just oscillating periodically at a frequency uj x o ~ ( e a;,o) ■ So, as far as the 
uncoupled dynamics is concerned, a typical state of a chain of identical anharmonic oscillators on the one 
hand, and a typical state of a disordered chain of harmonic oscillators on the other hand, are qualitatively 
completely similar. In the first case however, typical means typical with respect to the Gibbs measure, 
while in the second it means typical with respect to both the Gibbs measure and the disorder. 

However, when trying to adapt our strategy, the problem of resonances cannot be ruled out so easily. 
Indeed, in the disordered chain, resonances only affected some fixed sites, while here, the places where 
resonances occur also move with time, a situation which in fact favors the transport of energy. Still 
our main intuition remains that, due to quick averaging of fast oscillations, energy gets trapped in finite 
portions of space for time-scales that might not be an inverse polynomial of e. 
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